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Abstract. Calculation of gravitational radiation during 
binary inspiral leading to possible formation of a Thorne- 
Zytkow (TZ) object (a neutron star inside a supergiant 
core) is performed. The calculations were done for poly- 
tropic density distributions with different indeces n, as 
well as for realistic models of supergiants. A maximum 
frequency of the emitted gravitational waves during this 
process was found to range from a few to 300 initial kep- 
lerian orbital frequensies, that is from 10 -5 to ~ 0.1 Hz, 
depending on the model. A dimensionless strain metric 
amplitude can reach h ~ 10~ 235 for a source lying 10 kpc 
away from the Sun. We conclude that TZ objects forming 
at a rate of 1 per 500 yrs in the Galaxy could be poten- 
tial astrophysical targets for space laser intereferometers. 
Analysis of gravitational waveforms emitted during TZ- 
object formation could bring a unique information about 
stellar structure. 
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An important result of the extensive numerical calcu- 
lations is that the ratio of binding energy of the stellar 
envelope ejected during the CE stage to change of the 
system's orbital energy, o>ce = AE env / AE or i, is of or- 
der unity. This implies a more than 100-fold shrinkage of 
the initial orbital separation after the spiral-in for typical 
parameters of high- mass X-ray binaries (HMXB) (a super- 
giant with a mass M sg ~ 10 — 15Af Q having a ~ 2 — 3M Q 
helium core in pair with a m — 1.4M neutron star). 
Thus, a crude lower limit for formation rate of TZ-objects 
can be taken as that for HMXB, i.e.~ 10 -3 per year per 
a 1O 11 M Milky- Way-type galaxy. In fact, the formation 
rate of TZ objects in the Galaxy typically can be nearly 
twice as high. This means that > 30 formations of TZ ob- 
jects per year could be observed from within a distance of 
30 Mpc. 

Let us estimate the frequency and amplitude of gravi- 
tational radiation emitted during this process. Consider a 
close binary system consisting of a neutron star orbiting 
a massive red (super)giant. The dynamical friction is the 
main drag force acting the neutron star inside the super- 
giant's atmosphere: 



1. Introduction 

According to the modern scenario of binary star evolution 
(see van den Heuvel (1994) for a review), massive binary 
systems with mass ratios far from unity (>~ 2.5) or in 
wich the mass losing star has a convective envelope must 
pass through a so called common envelope (CE) stage 
(Paczyhski 1976). Considered for the first time in detail 
by Taam, Bodenheimer and Ostriker (1978), the the dou- 
ble core hydrodynamical evolution in common envelope is 
repeatedly calculated (see Livio 1994 for a review). The 
issue of the common envelope evolution is either ejection 
of part or all of the envelope, or coalescence of the two 
cores with little mass ejection leading to formation of a 
Thorne-Zytkow object (Thorne and Zytkow 1977). 



F d (r) ^ R 2 G p{r)v r (rf , (1) 



where r is radial distance to the supergiant's center, p is 
its density, v r is the neutron star velocity relative to the 
envelope, and Rg = 2Gm/ (v 2 + a 2 ) is gravitation capture 
radius for the neutron star (G is the newtonian gravity 
constant, a s is the sound velocity). 

Hydrodynamical calculations (Livio 1994) have shown 
that the neutron star can reach the dense stellar core dur- 
ing some 10 3 years and the final plunge inside the dense 
core takes less than one month. A crude estimation for 
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a dimensionless strain metric amplitude of gravitational 
radiation at the beginning of the plunge is 

r fl (M)r(m)^ / M \ 

RcD \10M Q J 

( m \ /10 u cto\ flOMpc\ ( ' 
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where M and R c are mass of the star and radius of its core, 
respectively, r g (M) = 2GM/c 2 is gravitational radius of 
mass M (c is the speed of light) and D is a distance to 
the object. Obviously, the amplitude must vanish as the 
neutron star approaches the center of the star. 

The maximum gravitational radiation output during 
the spiral-in is expected to occur at two-folded circular 
frequency of the neutron star orbital motion. Here we will 
calculate gravitational wave emission as the neutron star 
plunges into a model star with a polytropic density distri- 
bution. We will also consider the case of a realistic density 
profile in giant stars. 



Fig. 1. Evolution of gravitational wave frequency for different 
polytropic configurations. The top and right frames are nor- 
malized for a 2Mq core with R = 10 10 cm. 



2. Quadrupole gravitational radiation from 
spiralling-in neutron star 

The transverse-traceless components of the metric pertur- 
bations are connected with the second time derivatives 
of the transverse-traceless components of the system's 
quadrupole moment, (see Misner, Thorne, Wcelcr 

1973): 



h 
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so we need to calculate how the changes along the 
neutron star trajectory. 

2.1. Quadrupole moment of the system 

We consider a point mass m in a circular orbit with a 
radius r(t) < R inside a spherical star of mass M and 
radius R. For our model calculations we will assume the 
inner structure of the giant star not to change during the 
spiral-in. Lagrangian mass comprised inside the orbit is 
M(r) = 47T J r p(x)x 2 dx. 

The traceless part of the second moment of the sys- 
tem's mass density p computed in a Cartesian coordinate 
system centered on the binary barycenter is 



SSij(t) = / P(t){XiX 



S i jX l x l )d 3 x . 



(4) 



We assume the binary orbit to lie in the XY coordinate 
plane. We introduce a polar angle ip between the line con- 
necting centres of the bodies and OX axis. Symmetry rela- 
tive to the orbital plane implies Q xz = 3j, z = 0. Assuming 
the neutron star does not influences the matter distribu- 
tion inside the supergiant (which can be the case for close 
enough spiralling-in binaries, see Taam et al. (1978)) and 



the giant star rotation as a whole around the barycenter, 
the non-zero components of the quadrupole moment are 



%x = (cos Lp - -)na 



^s yy = (sin 2 (p-^-)pa 
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This is the well-known result for two point masses in the 
orbit. 

2.2. Trajectory of the neutron star 

Equations of motion of the neutron star (treated as a point 
mass) subjected to the dynamical friction force can be 
written as 



d_ f dL\ dL dD 



(10) 



Here L is Lagrange's function of the point in a gravity 
field 
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Fig. 2. Evolution of the gravitational wave frequency faw for Taam (1979) (a) and Weaver et al. (1979) (b) red giant models. 



and D is a dissipative function: 

D = R 2 G pv 3 = R 2 G p(r) [r 2 + r V] 3/2 



(12) 



Energy losses due to dynamical friction only are taken 
into account in the dissipative function. The gravitational 
radiation reaction contributes unsignificantly as it is of the 
order (v/c) 5 . We assume no change of the neutron star 
mass (see Chevalier (1993) for possible issues of accretion 
onto spiraling-in neutron star). 

The dynamical equations were numerically solved in 
order to obtain arrays of r(i), ip(t) and their time deriva- 
tives up to the second order, which were then used for 
computing the second time derivatives of the traceless 
quadrupole moment 

The density distribution inside the core was modelled 
by polytripic distributions with indeces n (so that the 
pressure is P oc p 1+1 /™) ranging from to 2.5. We also 
apply our calculations to a realistic density profile inside 
a 16 Mq red giant (Taam 1979), and to the extreme case 
of density distribution inside a 15 Mq supergiant just 
prior to the core collapse (Weaver et al. 1979), which ap- 
proximately corresponds to a polytropic distribution with 
n = 3. 

In Fig.[j] we present the evolution of the two-folded or- 
bital frequency (the main frequency of the gravitational 
wave radiation) for the neutron star, 



and the longer time the star remains in a keplerian orbit 
at the outer less denser layers. The oscillations seen for 
low n reflect elliptical character of the initially circular 
motion suddenly subjected to a drag force. Note that for 
polytropes with n > 1 the plunge to the final frequency 
occurs very rapidly, so the gravitational wave radiation 
will have a burst-like character for dense configurations 
(see the next section). 

Evolution of the gravitational wave frequency few 
with time for the case of density distribution inside Taam 
(1979) and Weaver et al. (1979) red giants is shown in 
Fig.§. 



(13) 



as a function of time (in units of the initial keplerian time) 
for different polytropes. Here tp denotes the dimcnsionless 
time derivative. The top and right frames are shown for 
M = 2M , R = 10 10 cm. It is clearly seen from the fig- 
ure that the larger n, the higher final frequency can be 



Fig. 3. Trajectory of the neutron star final plunge in a red 
giant with extreme density distribution (Weaver et al. 1979). 



The neutron star trajectory inside the innermost re- 
gion a < 0.2R for Weaver et al.'s red giant is shown in 
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Fig.||. The duration of the final plunge is of the order of a 
few 0.001 the initial keplerian time. 

2.3. Gravitational radiation 

The dimensionless waveforms h + and h x of the metric 
strain were calculated in newtonian quadrupole approx- 
imation (see e.g. Kochanek et al 1990). As an example, 
we calculated h+ and h x as seen from a distance D for 
the face-on orientation of the binary relative to the line of 
sight, so that 

Dh+ = ^ xx -% y ); (14) 



Dh> 



L 4 ^xy ■ 



(15) 



By measuring mass in units of the stellar mass M, radial 
distance in units of the stellar radius R, time in units of 
the initial free-fall time 



Fig. 4. Averaged over binary orientation metric strain ampli- 
tude (h 2 ) (in units of (r g /D) 2 ) from a neutron star spiraling-in 
inside a stellar core with a polytropic density distribution, as 
a function of time in units of v/ R 3 /GM. 



^.0^-1.6x103(^1^, 
we can rewrite these equations in the form 



r g (M) 



ARM 



D 



r g {M) 



l r g (M) M 
2 R M 



y i 



(16) 



(17) 



where ^>ij denotes the dimensionless components of the 
quadrupole moment Ssjj (cf. intuitive estimation @). 

In addition, we calculated averaged over orientation of 
the binary waveforms (Kochanek et al. 1990) 
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Temporal behaviour of these quantities for polytropes 
are shown in Fig.|| For more concentrated configurations 
(n > 1), a burst-like behaviour corresponding to the rapid 
final plunge is clearly seen. 

Fig.p and ^ shows the h + polarization of the gravi- 
tational wave calculated for a face-on orientation of the 
system relative to the line of sight, and the corresponding 
angle averaged metric strains (h 2 ) for cases of Taam's and 
Weaver et al's red giants. 



2.4- A more realistic treatment of the neutron star spiral- 
in 

We have assumed so far that the stellar structure is not 
changed by the spiralling-in neutron star, which is of cause 
a very crude zero approximation used in order to get sim- 
ple estimations of the effect. However, the qualitative be- 
haviour of the frequency and amplitude of generated grav- 
itational radiation is expected to not differ significantly 
in more accurate calculations. The problem itself is of 
course a complicated three-dimensional hydrodynamical 
one, but as an illustration we performed more precise 1.5- 
dimcnsional calculations of the neutron star spiral-in in 
case of Taam's model. 

We assumed after Taam (1979) that the entire com- 
mon envelope is rotating rigidly (by action of effective 
convective angular momentum transfer) with a spatially 
constant angular velocity defined at each moment as v e = 
JefNs/Ie, where J e is the total angular momentum trans- 
ferred from the orbit to the envelope rotation, v e the en- 
velope circular velocity, r/vs the position of the neutron 
star and I e the moment of inertia of the common envelope. 
The relative velocity of the neutron star is v r = v^s — v e- 



Thus, instead of (12) the dissipative function is 



D 



p{r)v 3 r 



(20) 



The precise value of sound velocity a s is a compli- 
cated problem and can be accurately calculated only under 
fully hydrodynamical treatment. We assume after Cheva- 
lier (1993) the motion of the neutron star in the envelope 
to be mildly supersonic with a Mach number of the or- 
der 3 — 1.4 in the central regions. Therefore, during the 
calculations we neglect by the sound velocity in the outer 
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Fig. 5. h+ - polarization of the gravitational wave (a) and the corresponding averaged over binary orientation metric strain 
amplitude (h 2 ) (b), emitted during spiral-in of a 1.4 M Q neutron star inside Taam's (1979) model of red giant density distribution 
(M = 15M Q , R = 4.2 x 10 12 cm), as a function of time in units of ^/R?/GM. 



Fig. 6. The same as in Fig.^j for Weaver et al.'s (1979) density distribution inside a 16M© red giant with a radius of 3.9 x 10 1 



layers of the giant, and take it to be a s ~ v r in the inner 
regions after v r has been equal to the keplerian velocity 
for the Lagrangian mass M(rjvs). 

As expected, the calculations yielded similar depen- 
dences of few an d h on time (see Fig .0a, b), with the 
major difference being in time the neutron star takes to 
reach the final plunge. In this case that time is much larger 
(as the dragging force is much less), and is about 10 3 years, 
in agreement with the cited hydrodynamical calculations. 
The final plunge starts at a distance from the center of 
15% the initial radius, and lasts about several days. The 
GW frequency and amplitude prove to be higher than in 
the simplest model case considered above by a half an 
order and one order of magnitude, respectively. 



3. Discussion and conclusions 

In the present paper we calculated gravitational wave 
emission during formation of a Thorne-Zytkow object - 
a red giant star with a neutron star core, which must be 
formed under certain conditions during evolution of mas- 
sive binary systems. 

To estimate the gravitational wave emission, we made 
a number of simplifying assumptions: 

1. no (or weak) influence of the spiralling- in neutron star 
to the density distribution inside the giant star; 

2. the dynamical friction is the only dragging force; 

3. no accretion induced mass change of the neutron star. 

The first assumption is the most severe, as during the 
spiral-in phase the whole envelope of the red giant can 
be lost. However, for enough tight initial configurations 
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Fig. 7. The frequency (a) and averaged over binary orientation metric strain amplitude (h 2 ) (b) of GW emitted during spiral-in 
of a neutron star in the red giant model by Taam (1979) calculated using more realistic assumptions described in the text 



(case of Taam's model) the density profile inside the star 
changes no more than an order of magnitude to the end 
of the spiral-in stage (Taam 1979). This means that under 
more realistic conditions the neutron star takes more time 
to reach the red giant's centre, so our estimates should 
be considered as lower limits. This is illustrated by calcu- 
lations of the neutron star spiral-in inside Taam's (1979) 
red giant model assuming rigid rotation of the common 
envelope. 

The second assumption seems to be justified by exten- 
sive numerical calculations of the common envelope evo- 
lution. The third assumption is valid while the radiation 
pressure prevents matter from accreting onto the neutron 
star at a rate exceeding the Eddington limit ~ 10 _8 M Q /yr 
(see discussion of another possibility in Chevalier 1993). 

As is seen from Fig.|[ the maximum value of the di- 
mensionless strain metric amplitude h during neutron star 
spiral- in toward the center of a compact (~ Rq) poly- 
tropic configuration can reach ~ 10~ 24 at a frequency of 
~ 10 ~ 2 — 10 _1 Hz, assuming the source lying at 10 Mpc 
away from the Sun. Real density distributions yield ap- 
proximately the same amplitudes but at much lower fre- 
quencies ~5x 1(T 5 - 2 x 1(T 4 Hz (Fig.|,0), which can 
only marginally enters into the sensitivity range for the 
planned LISA detector (Laser Interferometer Space An- 
tenna; see Hills & Bender (1995) for more detail). 

However, the situation is not excluded when the neu- 
tron star enters a dilute envelope of a 16 M Q envelope, 
makes it to be dynamically unbound but does penetrate 
into the inner 2-3 Mq dense core to form a TZ object. 
Then the frequency of the GW could well fall into the 
range 10~ 3 — 1CP 2 , where the LISA detector is the most 
sensitive. 

To conclude, we stress the importance of hydrodynam- 
ical modeling of processes considered in the present paper. 
Analysis of gravitational wave signal emitted during for- 



mation of a TZ-object can bring a unique direct informa- 
tion about density distribution inside the star, which is 
not available by other methods. 
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